Top jets as a probe of degenerate stop-NLSP LSP scenario in the framework of 

cMSSM 

Kirtiman Ghosh^, Katri Huitu^, Jari Laamanen'^ and Lasse Leinonen^ 
^'^''^ Department of Physics, and Helsinki Institute of Physics, FIN-OOOI4 University of Helsinki, Finland and 
^ Theoretical High Energy Physics, IMAPP, Faculty of Science, 
Radboud University Nijmegen Mailbox 79, P.O. Box 9010, NL-6500 GL Nijmegen, The Netherlands 
^ kirti.gh@gmail. com, ^ katri. huitu@hel.sinki.fi, ^j. laamanen@science.ru. nl, ^ lasse. lemonen&helsinki.fi 

The degenerate stop NLSP and neutralino LSP scenario is well motivated but hard to detect 
at the collider experiments. We propose a novel signature for detection this scenario at the Large 
Hadron Collider (LHC) and demonstrate its feasibility. 



Supersymmetry (SUSY) is one of the most attractive 
extensions beyond the standard model (SM) and search 
for SUSY is one of the prime goals of the Large Hadron 
Collider (LHC). The production of colored SUSY par- 
ticles, such as squarks and gluinos, at the LHC plays 
the most important role in discovering SUSY. SUSY sig- 
nature should show up best in multi-jets plus missing 
transverse energy {I^t) channel in i?-parity conserving 
scenarios. Ifx results from the weakly interacting sta- 
ble lightest SUSY particles (LSPs) appearing in the final 
states. With 4.7 (4.98) fb^^ integrated luminosity {^Cdt) 
at 7 TeV center-of-mass energy, no SUSY signature has 
been detected in the multi-jets plus Sfr channel by the 
ATLAS 4,] (CMS (aj) cohaboration. As a result, limits 
are imposed on the squark and gluino masses in the con- 
text of constrained minimal supersymmetric extension 
of the Standard Model (cMSSM) fi|. The phenomenol- 
ogy of cMSSM is determined in terms of five parame- 
ters, namely, common scalar mo and gaugino rni/2 mass, 
tri- linear mass parameter Aq, ratio of Higgs vacuum ex- 
pectation values tan/3 and sign of Higgs mixing param- 
eter sign(/i). In the context of cMSSM with tan/3 = 10, 
Aq = and /i > 0, for large toq, gluino mass below about 
900 GeV is excluded at the 95% confidence level (C.L.) 
from ATLAS and CMS multi-jets plus I^t search. Equal 
mass squarks and gluinos are excluded below 1400 GeV. 

The main motivation for postulating the existence of 
SUSY is to stabilize the electroweak hierarchy against ra- 
diative corrections. In view of this motivation, the lower 
bounds on the squarks and gluino masses seems already 
rather large. However, in the radiative corrections to 
Higgs boson mass parameters, third generation (s)quarks 
are most important at one loop order. Moreover, the AT- 
LAS and CMS searches are not sensitive to the direct pair 
production of only third generation squarks Q. There- 
fore, stop masses of a few hundred GeV are still allowed 
and, in fact, favored by finetuning arguments. 

The stop can be the next-to-lightest supersymmetry 
particle (NLSP). This can be realized in the context of 
cMSSM because of large third-generation Yukawa cou- 
pling. The off-diagonal trilinear term is proportional to 
Yukawa couplings. Therefore, it can generate large mix- 
ing and hence, large mass splitting in the stop sector and 
consequently produce a light stop (ti) as NLSP. 



There are several phenomenological motivations to be 
interested in the ti-NLSP. The ti-NLSP scenario is well- 
motivated to explain the relic density (RD) measured by 
WMAP d]. In the cMSSM, the lightest neutrahno (x?) 
is mostly binolike and thus, gives rise to RD which is 
too large. The observed relic density can be explained 
if the ii-NLSP coannihilate with the binolike Xi-LSP. 
Furthermore, the electroweak (EW) baryogenesis favors 
a light-fi In fact, a SUSY scenario with a binolike 
X^-LSP and a fi-NLSP can simultaneously explain the 
observed baryon asymmetry in the universe as well as 
the RD. Finally, cMSSM with fi-NLSP is also favored in 
view of a Higgs boson mass {mh) about 125 GeV observed 
by the LHC [7]. In order to get a Higgs boson around 
125 GeV, we must have either significant stop mixing or 
a large stop mass scale. Large stop mixing in cMSSM 
with ^i-NLSP may give rise to mh around 125 GeV. 

Collider searches for ti-NLSP at the LHC is straight- 



forward as long as 



is sufficiently large. In this 



case, ti signature should pop up in multi-jets plus I^t 
search. If the fi and Xi mass splitting is small (which is 
favored by RD) then the decay channels ii — >• hWx\ and 
tXi are kinematically forbidden and the four-body decay 
ii ff'bxi is strongly phase space suppressed. The 
dominant decay mode will be loop induced two-body de- 
cay ii — )■ cxi- However, due to small ii-Xi mass splitting, 
the charm jets as well as I^t will become very soft at the 
collider. Therefore, di-jet -I- I^t signature resulting from 
titi production will be buried deeply in the background 
events and cannot been selected out with the current 
LHC search approaches [1] . Several alternative solutions 
to detect signatures of light ii nearly degenerate with the 
X? have been proposed, ti pair production in association 
with two 6-jets single photon [13] or jet 11 1 have been 
studied. Stop pair production with a hard jet is the most 
promising channel and can probe m^^ upto 250 GeV [11 1 
at the LHC with ^/s = 14 TeV and JCdt =100 fb~^ 

In this letter, we propose a novel collider signature 
for light stop search. In cMSSM with fi-NLSP, gluino 
(g) is heavier than ?i by a factor ~ 5 and thus, at the 
LHC, gg production cross-section is atleast two orders of 
magnitude smaller than the titi production cross-section. 
However, we find that gg production gives rise to a in- 
teresting signature which could be observed with a few 



2 



fb~^ luminosity of the LHC running at -y/s = 14 TeV. 

Gluino is lighter than the first two generation squarks 
(q) in most part of the cMSSM parameter space. If ti is 
the NLSP then g will dominantly decay to the tti pairs 
followed by ti — )■ cxi decay. In this case, gg produc- 
tion gives rise to two top quarks, two charm quarks and 
l^T- PP ^ gg ^ tt tltl — >■ tt cc XiXi- Since nig ^ m^^, 
the top and the stop arising from the gluino decay will 
be highly boosted. However, the charm quarks resulting 
from the decay of the boosted stops will be soft because of 
their small mass and small mass difference between stops 
and LSPs. The momentum of the boosted stops will be 
carried away by the massive lightest neutralinos. Top 
quarks further decay hadronically {t — bqq') or lepton- 
ically {t ^ blv). However, the highly boosted tops will 
decay to collimated collections of particles that look like 
single jets. It has been already established that top jets 
can now be considered as standard objects for event anal- 
ysis at the LHC [13,113. In this letter, we have studied 
two high transverse momentum (px) top jets plus large 
as a signature of cMSSM with fi-NLSP. 

Gluino pair production in this scenario also gives rise 
to usual multi-jet -I- i^T (when tops decay hadronically), 
opposite-sign dilepton-|-i^r(when gg decays to tt and tops 
decay leptonically) and same-sign dilepton (SSD)-l-i^T 
[l3 | (when gg decays to same sign tops and tops decay 
leptonically) signature. However, all these signatures suf- 
fer from the SM tt and QCD multi-jet backgrounds. Two 
high pt top-tagged jets-|- i^T signature has the following 
advantages: (i) The mis-tagging efhciency of gluon or light 
quark jets to be tagged as top jets is very small [13| . 
Therefore, top tagging significantly reduces SM multi-jet 
background, (ii) Large pt cuts could be applied on the 
top-tagged jets to reduce the SM tt background, (in) Fi- 
nally, this is a characteristic signature of cMSSM with 
ii-NLSP which is well motivated from WMAP data, EW 
baryogenesis and the LHC indication of ruh ~ 125 GeV. 

In Ref. [l^, ATLAS collaboration has presented a 
search for new phenomena in tt + I^t events at the LHC 
with ^/s = 7 TeV and JCdt =1.04 fb~^. In view of consis- 
tency of the data and SM expectations, 95% C.L. upper 
limit of 1 . 1 pb is set on cross-section times branching ra- 
tio of any new physics process which gives rise to tt + I^t 
signature. Similar signature results from gg production 
in cMSSM with fi-NLSP. Assuming 100% branching frac- 
tion for g — > tt*, naive estimation of next-to-leading or- 
der (NLO) gg production cross-section implies an lower 
bound about rrig ~ 560 GeV in this scenario. 

In cMSSM, ti can be the NLSP when the mixing term 
of the left and right stop states, = nitiAt — ^cot/3), 
is large. The largest mixing occurs when fi and At have 
opposite signs. In our analysis, we consider > and 
thus, scan over negative values of Aq. In addition to the 
mixing, the renormalization group running plays an im- 
portant role in the determination of the NLSP. In the 
scalar RGEs, the terms proportional to the Yukawa cou- 
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FIG. 1. The 3-dimensional parameter space gives rise to stop 
NLSP, for tan/3 = 15 and ^ > 0. To select a particular value 
for ^0, we demand m^^ — m-o ~ 35 GeV. The white dots 
corresponds to the points resulting into a ruh ~ 125 ± 2 GeV. 



TABLE I. Benchmark parameters in GeV (tan/3 — 15, fi > 0). 



BPl 


BP2 


BP3 


mo 


"^l/2 


Ao 


mo 


'"1/2 




mo 






1700 


420 


-3763 


1800 


450 


-4001 


1950 


491 


-4349 



plings (and scalar masses) decrease the soft masses, while 
the terms proportional to the gauge couplings (and gaug- 
ino masses) have an opposite effect. Due to the largeness 
of the strong coupling constant, squark masses tend to 
increase more than the slepton masses, even though the 
third-generation squarks have large Yukawa couplings. 
Thus, a small toi/2 compared to toq is preferred in or- 
der to suppress the strong coupling term in the squark 
RGE, which then may lead to a stop NLSP instead of 
stau NLSP. In our parameter space scan, we have varied 
Too over 1 to 2 TeV and TOi/2 over 200 to 800 GeV. Fig.[T] 
shows a 3-dimensional scan over toq (along x-axis), mi/2 
(along y-axis) and Aq (color gradient) for tan/? = 15 and 
fi > 0. In this parameter space, ti is always NLSP. For a 
given mo and mi/2, there is a range of ^o values which 
give rise to ti-NLSP. In Fig. [l] we consider the partic- 
ular ^0 (motivated form the observed WMAP RD [li|) 
35 GeV. The particle spectrum 



for which 



TO^o 



was calculated using SuSpect [17| with nit = 173.2 GeV, 
TOfc = 4.2 GeV, mr = 1.777 GeV and Usimz) = 0.1172. 

For this parameter space, we have checked the follow- 
ing constraints: For b ^ sj branching ratio, we consider: 
BRiB Xs-f) = (355 ± 142) x 10-^ including the the- 
oretical uncertainty in the SM as well as in the MSSM 
[13 ■ We have also checked that Bg — > fi^fi'^ branching 
ratio is smaller than the experimental upper limit [lij : 
BR{Bs ->■ < 4.3 X 10-8 (95% C.L.). The new 

physics (NP) contribution to the BR{B~^ — > r+z/^) can 

be quantified by defining: R^J^ = ""SB+'-Ir+^O^r" ' 
The 95% C.L. allowed range for R^^ is 0.99 < R^f^ < 
3.19 |20|]. With the combined data from WMAP, BAG 
(Baryon Acoustic Oscillations) in the distribution of 
galaxies and observation of Hubble constant, the density 
of cold dark matter in the universe is determined to be [5| 
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TABLE II. Masses (in GeV) of relevant sparticles. 





g 


UL 


UR 


t2 


ii 


b2 


bi 


xf 


X? 


BPl 


1064 


1880 


1872 


1334 


217 


1776 


1287 


352 


181 


BP2 


1133 


1994 


1985 


1411 


229 


1883 


1365 


378 


194 


BPS 


1227 


2161 


2150 


1523 


250 


2040 


1478 


414 


213 



TABLE IV. Signal and SM background cross-sections at the 
LHC with a/s = 14 TeV after different set of cuts. 



TABLE III. Model contribution to different observables. 



BPs 


mu 




BR(B ^ X,7) 


BR(Bs m^M") 


n>NP 


BPl 


124.0 


0.115 


2.77 X 10"* 


3.14 X 10"^ 


0.997 


BP2 


124.1 


0.108 


2.82 X 10~* 


3.13 X 10"® 


0.998 


BP3 


125.3 


0.131 


2.89 X 10"* 


3.12 X 10"® 


0.998 



n^h^ = 0.1126 ± 0.0036. We consider the preferred RD 
range of 0.0941 < fi^/i^ < 0.1311 at 2a level. We have 
used micrOmegas (v.2.4.R) [21] for calculating RD and 
other constraints. Fig. [T] shows that large part of the pa- 
rameter space is below the WMAP upper limit (and can 
be made consistent with the preferred RD range) as well 
as consistent with other experimental observations. 

The ATLAS and CMS collaborations have pre- 
sented the first indication for a Higgs boson with m/j ^ 
125 ± 2 GeV. Due to the large stop mixing, the fi-NLSP 
scenario can accommodate such a Higgs boson mass. In 
Fig. [1] the white dots correspond to the points resulting 
into a rrih ~ 125 ± 2 GeV. Fig. [T] shows that small part 
of mo-TOi/2 plane is consistent with both ^ 125 ± 2 
GeV and WMAP data as well as other constraints. With 
this motivation, we have chosen three benchmark points 
(BP), listed in Table U for presenting our numerical re- 
sults. The masses of relevant sparticles for the BPs are 
presented in Table [III In Table IIIIl we present the Higgs 
boson mass, RD and other constraints for the three BPs. 

We consider boosted top jets as the signature of 
cMSSM with ii-NLSP. To detect top jets, we have used 
Johns Hopkins top tagger algorithm [l3| . A brief descrip- 
tion of the algorithm is as follows: Particles are clustered 
into jets of size R using the Cambridge- Aachen (CA) al- 
gorithm [i^l. To find subjets, each jet in the event is 
declustered into two objects by reversing the CA cluster- 
ing. The softer object is thrown out if the ratio of its pr 
and the full jet pt is less than a parameter 5p, and the 
declustering continues on the harder object. The declus- 
tering stops and the original jet is considered irreducible 
when the px ratios of both objects are smaller than i5p, 
two objects are too close |A?7| — |A^| < 5r or there is 
only one cell left. The original jet declusters into two 
subjets if the pr ratios of both objects is larger than 5p 
and the previous steps are repeated on those subjets. To 
reconstruct boosted hadronic top jets, cases with 3 or 4 
subjets are selected for the further analysis. 

In our analysis, we consider the events with total 
hadronic energy deposits, Et > 1 TeV and use FAST- 
JET v.3.0.2 [231] for the clustering and subjet analy- 
sis. Since more highly boosted tops will be more col- 
limated, following Ref. [13|, we choose R and 5 parame- 





SM Background cross-section in fb 


Process 


Acc. 


Acc. Cuts + 


Acc. Cuts -I- 


Acc. Cuts + 




Cuts 


^T> 100 GeV 


^T> 300 GeV 


^T> 500 GeV 


tt 


60.09 


1.29 


8.1x10"^ 


>1.5xl0"^ 


jets 


51.32 


0.26 


6.2x10"-^ 


>3.1xl0"^ 


tiZ 


0.0022 


1.6x10"^ 


3.2x10"* 


>3.4xl0"* 


ttW 


0.074 


2.7x10"^ 


2.8x10"^ 


9.2x10"* 


BPs 


Signal cross-section in fb 


BPl 


2.49 


2.32 


1.51 


0.76 


BP2 


1.75 


1.66 


1.19 


0.73 


BP3 


1.22 


1.16 


0.87 


0.57 



ters as: for Et > 1, 1.6, 2.6 TeV, R = 0.8, 0.6, 0.4 and 
Sp(^r) = 0.10(0.19), 0.05(0.19), 0.05(0.19), respectively. 
A jet is considered as top-tagged if (i) it contains 3 or 
4 substructure, (ii) the invariant mass of the sum of the 
subjets momentum four-vectors falls in the range 150-200 
GeV (~ nit ± 25 GeV) and (iii) there exist two subjets 
which reconstruct the W-mass to within 65-95 GeV. 

We have used Prospino 2.1 with CTEQ6.6M 
parton distribution function (pdf )for computing NLO gg, 
95; 99 and qq* production cross-sections at the LHC 
with -y/s = 14 TeV. The QCD factorization and renor- 
malization scales are kept fixed at the sum of masses of 
the produced SUSY particles. Signal events, showering, 
hadronization, initial state and final state radiation are 
simulated by PYTHIA 6.421 An event is selected 
for further analysis if it contains at least two top-tagged 
jets with Pt > 500 GeV and \r]\ < 2.5. We reject events 
with one or more isolated electron or muon [27| with 
Pt > 20 GeV and ]ry] < 2.5. We collectively refer this set 
of cuts as Acceptance Cuts. The signal is characterized by 
large I^t resulting from the boosted Xi which escape the 
detector without being detected. Signal cross-sections for 
the BPs are presented in Table IIVI after acceptance and 
different cuts (^t> 0, 100, 300 and 500 GeV). 

The dominant SM background arises from the tt pro- 
duction. The efficiency for gluon or quark jets to be 
mistagged as top jets is very small [l3|- However, due 
to huge production cross-section, QCD multi-jet events 
might also contribute to the background. In these contri- 
butions, i^Tarises from momentum mis-measurement due 
to finite detector resolution. We have approximated the 
detector resolution effects by smearing the jet energies 
with Gaussian functions [i^ . Table IIVI shows that large 
i^T cuts reduce ti and multi-jets backgrounds to a negligi- 
ble level. Production of tt pairs with a Z or Ty-boson may 
contribute to the background when Z decays invisibly 
and W decays leptonically. Lepton veto in the acceptance 
cuts significantly reduce the ttW contribution. Multi-jets 
and tt were simulated with PYTHIA. Whereas, ttZ and 
ttW events were generated using ALP GEN [29], and the 
unweighted event samples were analyzed by PYTHIA. 
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FIG. 2. Constant signal cross-section contours in mo-mi/2 
plane after Acc. cuts + ^t> 500 GeV at the 14 TeV LHC. 
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FIG. 3. Signal cross-sections as a function of stop mass for 
two different set of mo-mi/2 after Acc. cuts + i^T> 500 GeV 
at the 14 TeV LHC. tan/3 = 15 and /i > are assumed. 

For background analysis, we have used CTEQ6L1 [s^ 
pdf. The renormalization and factorization scales have 
been set equal to the subprocess center-of-mass energy 
For tt production, we have used NLO-f NNLL cross- 
section of 877 pb at 14 TeV LHC. In Table|lVl we present 
SM background contributions after different I^t cuts. 

According to Table HVl for I^t> 500 GeV, the number 
of predicted SM background event with JCdt = 100 fb^^ 
is less than one. If we demand 5(10) signal events for 
the discovery, BPl, BP2 and BPS will be probed with 
7(13), 7(14) and 9(18) fb"! integrated luminosity of the 
LHC running at -y/s = 14 TeV. In Fig. [2l we present con- 
stant signal cross-section contours on mo-TOi/2 plane af- 
ter Acc. cuts + ^T> 500 GeV at the LHC with = 14 
TeV. Assuming 5(10) signal events for the discovery, the 
0.7, 0.3 and 0.15 fb cross-section contours (equivalently, 
for large mo, nig ^ 1.1, 1.5 and 1.7 TeV ) can be probed 
with JCdt = 7(13), 17(33) and 33(67) fb^i. 

Up to this point, we have discussed two high-p^ top- 
tagged jets+i^T-as the signature of nearly degenerate ti- 
NLSP LSP scenario. To show the feasibility of this sig- 
nal for non-degenerate ti-Xi scenario, we have varied Aq 
and hence, mj-^, keeping mg and toi/2 fixed. In Fig. [31 
we present the signal cross-sections (for 14 TeV LHC) 
after Acc. cuts + I^t > 500 GeV as a function of m^^ . 



Initially signal cross-section decreases because increas- 
ing rrij^ implies less g-ti mass splitting and hence, less 
boosted tops. In Fig. [31 mi/2 = 700 GeV corresponds to 
m^o ^ 300 GeV. For TOj^ > 475 GeV, boosted ti decay to 
txi increases the possibility of getting more boosted top 
jets and thus, enhances the signal cross-section. Fig. [3] 
shows that signal cross-section does not vary much as 
long as stop is the NLSP. Therefore, the signature under 
consideration is equally applicable for the whole ti-NLSP 
region of the cMSSM parameter space. 

To summarize, we have investigated two high-pr top- 
tagged jcts+ij/T as the signature of cMSSM with degen- 
erate fi-NLSP LSP. We have considered gg production. 
Boosted tops results from g — > tti decay. We find that as 
long as stop is the NLSP (including degenerate ii-NLSP 
LSP scenario), collider searches for high-pT top jets plus 
large I^t can probe gluino mass upto 1.7 TeV at the 14 
TeV LHC with 33 fb~^ integrated luminosity. 
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